Magnetic and Transport Properties of Fe-Ag granular multilayers 
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Results of magnetization, magnetotransport and Mossbauer spectroscopy measurements of se- 
quentially evaporated Fe-Ag granular composites are presented. The strong magnetic scattering of 
the conduction electrons is reflected in the sublinear temperature dependence of the resistance and in 
the large negative magnetoresistance. The simultaneous analysis of the magnetic properties and the 
transport behavior suggests a bimodal grain size distribution. A detailed quantitative description 
of the unusual features observed in the transport properties is given. 
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I. INTRODUCTION 

As promising candidates for magnetic recording and 
sensor applications heterogeneous magnetic materials, 
i. e. multilayer structures^* 2 " of alternating ferromagnetic 
and nonmagnetic layers and granular composites*™*^ 
have been studied widely in the last two decades. The 
giant magnetoresistance (GMR) in these systems have 
been explained by elastic scattering of the conduction 
electrons on magnetic moments of differently aligned 
magnetic entities. Gittleman et al. have shown that 
in superparamagnetic granular alloys this considera- 
tion leads to a magnetoresistance proportional to the 
square of the magnetization** Deviations from this rela- 
tion has been attributed to the size distribution of the 
magnetic scatterersA*^ and interactions between these 
scatterer o 10 : 11 : 12 Various assumptions on the form of the 
size distribution have been made in order to obtain a 
phenomenological description of the GMR phenomena in 
different granular systems. 

In this paper we present a systematical study of 
the magnetic and magnetotransport properties of vac- 
uum evaporated granular Fe-Ag structures. The ob- 
served large, negative non-saturating magnetic field de- 
pendence and the unusual sublinear temperature depen- 
dence (d 2 R/dT 2 < 0) of the resistivity have been an- 
alyzed simultaneously. This allowed the separation of 
the various scattering processes and the identification of 
two characteristic grain size determining the macroscopic 
magnetic and transport properties without making any 
assumption on the grain size distribution. 



II. EXPERIMENTAL 

The Fe-Ag multilayer samples were prepared by 
sequential vacuum evaporation in a base pressure of 
10~ 7 Pa onto Si(l 1 1) single crystal substrates at room 
temperature. The mass of the deposited material was 
measured by a quartz oscillator and the nominal layer 
thickness was calculated using the bulk density of Fe and 
Ag. In this paper specimens prepared with the following 



sequences are discussed: 

(A) [Ag(2.6nm) /Fe(0.2nm)] 75 /Ag(2.6nm) 

(B) [Ag(1.3nm)/Fe(0.2nm)] 75 /Ag(1.3nm) 

(C) [Ag(0.8nm)/Fc(0.2nm)] 75 /Ag(0.8nm) 

Structural characterization of the samples by X-ray 
diffraction indicated a nanometer scale grain size of the 
constituents [for details see Ref. 13], however, due to 
the strong overlap between the diffraction lines of bec- 
Fe and fcc-Ag, a quantitative evaluation of the size of 
the magnetic Fe grains was not possible. The absence 
of peaks of the X-ray reflectivityii* when the Fe layer 
thickness is less than w 1 nm is attributed to discontinu- 
ities of the Fe layers. This limit of the continuous layer 
regime is quite similar to that observed in other multi- 
layer systems^*^**** composed of transition metals with 
immiscible nonmagnetic elements and these kind of dis- 
continuous multilayers are often referred to as granular 
multilayersit The relation between the size of the grains 
in granular multilayers and the nominal layer thickness 
is generally determined by a three-dimensional growth 
processA^and the average diameter of the grains can be 
much larger than the nominal layer thickness. It depends 
on the material parameters (e.g. lattice parameter mis- 
match, surface energy, etc.), as well as various parame- 
ters (deposition rate, substrate temperature, etc.) of the 
deposition technique. The decrease of the spacer layer 
thickness was foundA*^ to increase the magnetic grain size 
in an other series of Fe/Ag granular multilayer samples. 

The magnetic structure of the samples was exam- 
ined by a superconducting quantum interference device 
(SQUID) and by transmission Mossbauer spectroscopy. 
For the latter purpose the samples have been removed 
from the substrate and folded up to achieve an appropri- 
ate thickness for transmission measurements. 

The largest GMR and the strongest magnetic scatter- 
ing of the conduction electrons were observed in sample 
A, as it will be shown later. The 57 Fe Mossbauer spectra 
of this sample taken at several temperatures are shown 
in Fig. ^ The room temperature spectrum contains a 
paramagnetic doublet with a large isomer shift relative 
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FIG. 1: Transmission Mossbauer spectra of sample A at vari- 
ous temperatures. At T — 4.2 K spectra with and without an 
applied magnetic field are shown. 

to cv-Fe (0.18 mm/s) and a quadrupole splitting (0.45 
mm/s) characteristic to a system of small Fe clusters em- 
bedded in an Ag matrij*2S. The low temperature spectra 
show that the sample is superparamagnetic (SPM) and 
as the magnetic clusters gradually freeze below the block- 
ing temperature (around 50K), the six-line pattern char- 
acteristic to the Zeeman splitting of the nuclear levels 
of 57 Fe appears. At T = 4.2 K the paramagnetic frac- 
tion is absent, but the spectral lines are much broader 
and the hyperfine parameters are different than those of 
bulk bcc-Fe. The 4.2 K spectrum could be fitted with a 
distribution of hyperfine fieldsSi with an average value, 
H av — 34.7 T and standard deviation cth =2.4 T. 

It is worth noting that the intensity ratios of the six 
lines indicate a significant spontaneous alignment of the 
magnetic moments. The intensity of the six lines of a 
sextet is distributed as 3 : I2-5 : 1 : 1 : 72-5 : 3, where 
L2-5 = 4 sin 29/(1 + cos 20) is the intensity of the m = 
transitions, and O is the angle between the direction of 
the gamma-ray (perpendicular to the sample plane) and 
the magnetic moment. In case of a random distribu- 
tion of the magnetization directions I2-5 — 2. The ob- 
served small intensity, I2-5 — 0.5, indicates a close-to- 
perpendicular alignment of the magnetic moments with 
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FIG. 2: Zero field cooled (ZFC) and field cooled (FC) magne- 
tization of samples A, B and C as a function of temperature 
[panels (a), (b) and (c) respectively]. 



respect to the sample plane. Applying a magnetic field 
perpendicular to the sample plane could fully align the 
moments parallel to the field, as it is indicated by the 
I2— 5 = intensity. 

The spectra measured in external magnetic field were 
also fitted with a distribution of the hyperfine fields and 
the parameters obtained are: H av = 32.5 and 27.8 and 
<7h = 2.3 and 2.4 T for B cxt = 3 and 7T. The external 
magnetic field does not affect the width of the distri- 
bution indicating that it is not due to relaxation of the 
magnetic moments. The broad spectral lines of the ob- 
served sextet result from a distribution of the Fe neigh- 
borhoods, which can be due both to the large number 
of surface atoms in small grains and to non-equilibrium 
mixing of the elements^ during the growth process. The 
field-independent width and the decrease of the observed 
average hyperfine field in external field indicate the fer- 
romagnetic alignment of the magnetic grains. The hy- 
perfine field of a ferromagnet is decreased by the applied 
field, since it is oriented antiparallel to the magnetic mo- 
ment. We note that at 4.2 K the statistical errors allow 
an SPM fraction containing less than 2 atomic % of the 
Fe atoms. 

The freezing of the superparamagnetic moments - seen 
in the Mossbauer spectra - also appears in the tempera- 
ture dependence of the low field susceptibility measured 
by the SQUID. Figure |2l shows the results for the 3 sam- 
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pies, after cooling them either in zero or in 1 mT per- 
manent magnetic field. The blocking temperatures of 
about 40, 150 and 300 K for samples A, B and C, respec- 
tively can be deduced from this experiment. Sample A 
(Tb = 40 K) exhibits a textbook example for the super- 
paramagnetic behavior. The maximum of the ZFC curve 
is more smeared out for samples B and C and in case of 
sample B the irreversibility temperature, where the FC 
curve starts to deviate from the ZFC one, is much larger 
than the Tg defined by the ZFC maximum. Such fea- 
tures are generally explainedSiSi as the effect of a grain 
size distribution and interaction among the magnetic par- 
ticles, which obviously play a role as the average grain 
size and the concentration of the magnetic component 
increase. 

The magnetization versus applied magnetic field curves 
of sample A are shown in Fig.|31up to B = 5 T at four dif- 
ferent temperatures. Above the blocking temperature the 
magnetization can quite well be described^ by a single 
Langevin function, as expected for a superparamegnetic 
system of uniform grains. The average magnetic moment 
of the SPM grains deduced from the fit is = 500 ^b • 
It indicates that the bulk magnetic properties are mainly 
determined by the average size grains (about 1.8 nm) and 
the effect of the size distribution is negligible. 

The magnetotransport measurements have been per- 
formed in the current in plane arrangement by four con- 
tact method. The coincidence of the parallel (H || /) 
and transversal (H _L I) magnetoresistance and the ab- 
sence of any anisotropic component is characteristic to 
the GMR phenomenon in granular systems^. The tem- 
perature dependence of the resistivity of samples A — 
C in zero and 12 T magnetic field normalized to their 
T = 290 K values are displayed in Fig. Contrary to or- 
dinary metallic systems, the resistivity is sublinear above 
40 K for each sample. Similar observation was reported 



300 




2 3 
Magnetic field (T) 



FIG. 3: Magnetization of sample A up to B = 5T magnetic 
field at four different temperatures. The solid line on the 
room temperature curve confirmes the SPM Langevin model 
with a characteristic magnetic moment of 500 (ib ■ 
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FIG. 4: Resistivity of samples A-C (normalized to their T = 
4.2 K, B = 12 T values) as a function of temperature at zero 
(a) and 12 T magnetic field (b). 



by Milner et al. for granular systems^ but only in pres- 
ence of a high magnetic field. As discussed later in detail, 
the above qualitative feature of the zero field tempera- 
ture dependencies signifies the presence of an extremely 
strong magnetic scattering in our samples. 

The anomalous character of the p{T) curves is the most 
dominant in case of sample A. Simultaneously, this sam- 
ple exhibits the largest magnetoresistance at low tem- 
perature, as shown in the bottom panel of Fig. [S] The 
magnitude of the GMR measured at T = 4.2 K in a field 
of B = 12 T is 26%, 18% and 16%, for samples A, B 
and C, respectively. Note that at room temperature this 
order is reversed. 

The magnetoresistance curves of the three samples 
have a common feature, they do not saturate even in 
high magnetic fields and at low temperatures, where the 
magnetization seems already to be saturated. This in- 
dicates that significant magnetic scattering takes place 
at magnetic entities much smaller than the typical grain 
size which determines the macroscopic magnetization. It 
is well known that the scattering amplitude of ferromag- 
netic grains embedded in a nonmagnetic metallic ma- 
trix is size dependent and the contribution of the smaller 
clusters is strongly enhanced^ The magnetic moments of 
these smaller clusters are harder to rotate by an applied 
magnetic field thus the saturation of the magnetoresis- 
tance is slower than that of the net magnetization arising 
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FIG. 5: Magnetoresistivity of samples A-C up to 12 T mag- 
netic field at room temperature (a) and T = 4.2 K (b). 



dominantly from the larger grains. Note that similar be- 
havior was found found in Fe-Ag co-deposited granular 
filmsSi as well as in many other systems (for a review 
seeSi). In the next section we give a more detailed anal- 
ysis of the above qualitative picture. 



III. ANALYSIS 

In order to identify the resistivity contribution of 
the magnetic scattering process we assume that the 
Matthiessen-rule can be applied, i.e the resistivity is com- 
posed of 3 terms: 

p(T, B) = po + Pph(T) + p mag „(T, B) . (1) 

Here po denotes the residual resistivity, p p h(T) is the con- 
tribution arising from the phonon scattering, and the re- 
maining part is attributed to the magnetic scattering. 
The separation of the latter term requires further as- 
sumptions. 

Instead of making assumption on the shape of the grain 
size distribution function^*! we use physical considera- 
tions to get insight into the size distribution of the grains 
determining the magnetic scattering. It is based on the 
analysis of the temperature dependence of the resistivity 
measured in zero and high magnetic field, as shown in 
Figs. HO and [7| In these limits the large grains - which 



dominate the bulk magnetic properties - do not give tem- 
perature dependent contribution to the magnetic scatter- 
ing, as the system is either nonmagnetic (zero field limit, 
above the blocking temperature) or fully polarized fer- 
romagnetically (high field limit). We will see that the 
small clusters can also be characterized by a single aver- 
age size, and in high fields the temperature dependence 
of the magnetic scattering is determined solely by this 
characteristic size: the magnetic moment of these small 
clusters is the only fitting parameter for the calculated 
curves describing the shapes of p magn (T, B = 12T) shown 
in Fig. EI Finally, the consistency of our analysis will also 
be demonstrated by evaluating the magnetic field depen- 
dence of the resistivity in the whole magnetic field range 
(Fig. [!5J using the characteristic sizes determined inde- 
pendently from the magnetization data and the temper- 
ature dependence of the resistivity in the zero and high 
field limits. 

In a granular system the magnetic scattering depends 
on the correlation between the localized magnetic mo- 
ments of the grains, (/2i/If), while a conduction electron 
is scattered from /2j to /2f within its spin-diffusion length' . 
In zero magnetic field, well above the blocking tempera- 
ture the magnetic moments of the grains (including the 
large grains) are fully disordered. In these circumstances 
no temperature dependence is expected from the mag- 
netic scattering, at least until the spin diffusion length 
is large enough. We assume that at high temperatures 
the temperature dependence of the resistivity arise solely 
from the phonon contribution, 

/ T \3 ,e/T 2 d 

As the the phonon term is linear above the Debye tem- 
perature (8 w 210 K)&, the strength of phonon scatter- 
ing, oi, can be determined from the high temperature 
slope of the zero field resistivity curves. The calculated 
Pph(T) curves are shown in Fig. [fj] by dashed lines for the 
3 samples. The difference of the total resistivity and 
its phonon related part is attributed to the magnetic 
scattering, and po + p ma , gn (T, 0) is displayed in Fig. 
by dotted lines for each sample. As it was expected, 
the magnetic scattering is temperature independent at 
high temperatures and the lower the blocking tempera- 
ture the wider is the flat part of the separated magnetic 
contribution curve. It is also seen that the magnetic 
scattering gradually decreases as the blocking temper- 
ature is approached from above. Note, however, that at 
T = 4.2 K there is still a considerable contribution from 
it, i.e. p(T = 4.2K, B = 0T) ^ p(T = 4.2K, B = 12T) as 
it can be seen from Fig. 

Since the phonon term is magnetic field independent, 
the /9 p h(T) curves determined from the zero field tem- 
perature dependencies can be used to separate the mag- 
netic scattering contribution in the high field measure- 
ments. This is shown in Fig. [7| for the B = 12 T mea- 
surements. A comparison of the related panels of Fig. 
and Fig. [7| reveals that sample A has the strongest nega- 
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FIG. 6: Analysis of temperature dependence of the resistiv- 
ity of samples A-C [(a) -(c), respectively] at zero magnetic 
field. The solid lines are experimental data, dashed lines indi- 
cate the temperature-dependent resistivity contribution aris- 
ing from phonon scattering as calculated from the high tem- 
perature slope of p(T) and the assumed Debye-temperature. 
The dotted curves are the differences of the former two shifted 
upwards by the residual resistivity, and are attributed to the 
temperature-dependent magnetic scattering on small Fe clus- 
ters. 



FIG. 7: Separation of the temperature-dependent resistivity 
contributions related to different scattering mechanisms at 
B = 12 T magnetic field. The magnetic field independent 
phonon part is indicated by dashed line, the magnetic scat- 
tering related contribution is plotted by dotted line for each 
sample [A-C in panels (a) -(c), respectively]. 



characteristic moment S is described by the Brillouin- 
function: 



tive curvature of the resistivity and the biggest change in 
Pmagn(r, B — 12T). In this sample magnetic scattering 
at B = 12 T dominates over even the phonon term in a 
very broad temperature range. 

In order to describe the evaluated /9 ma gn(T, B = 12T) 
curves we assume that in the high field limit the magnetic 
scattering of the spin-polarized electrons is proportional 
to the spin disorder of the small clusters. The magnetic 
moments of the large grains are fully aligned by the ap- 
plied magnetic field, as it could be deduced from the 
Mossbauer spectra of Fig. ^ The spin disorder for a 



p masn (T,B) = a 2 (S-(S z )) 
= a 2 
1 



-coth^ 
2 2kT 



(3) 



Here, S and (S z ) are the total spin and its z-component 
of the scatterers, respectively. 

The fitted p m agn(T, B — 12T) curves are shown in 
Fig. [HI by solid lines. Apart from a normalization factor, 
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FIG. 8: Temperature dependence of the resistivity contribu- 
tion attributed to magnetic scattering on tiny Fe clusters in 
B — 12 T magnetic field. Symbols are experimental data after 
subtracting the phonon contribution and as shown in Fig. |7| 
solid curves are calculated from Eq. |3] with fitting parame- 
ters S — 16.6 fiB, 17.0 /xb and 12.5 /^b for samples A-C, 
respectively. 



the only fitting parameter is the magnetic moment char- 
acteristic to the small Fe clusters. The good agreement of 
the experimental and the calculated curves indicates that 
the size distribution of these clusters is negligible. The 
fitted values are in the same order of magnitude for all 
the samples; S — 16.6 /ie, 17 /ib and 12.5 /ie for samples 
A-C respectively. 

Next we discuss the magnetic field dependence of the 
resistivity of sample A, which exhibits the strongest mag- 
netic scattering in the superparamagnetic state, where 
the process of magnetic saturation is well understood. 
For a numerical analysis we use the two characteristic 
sizes determined from the previous experiments. The 
Langevin fit to the magnetization experiments performed 
in the superparamagnetic temperature range has shown 
the presence of large grains with ^g ~ 500 /^b (Fig. GJ), 
while the temperature dependence of the resistivity in 
high magnetic field indicated the presence of small clus- 
ters with S « 17 hb- Following Gittleman's modefi we 
describe the magnetoresistance by the field dependence of 
the correlation between the localized magnetic moments 
responsible for an initial and a final magnetic scattering 
process: 

— oc (flipf) oc (jkB) (pfB) ■ (4) 
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FIG. 9: Magnetoresistivity of sample A in the superparam- 
agnetic range, at room temperature. The open symbols are 
experimental data, the solid line is the fit according to Eq. 
The dashed and dash-dotted curves correspond to the first 
and second terms of Eq. respectively. The magnetic mo- 
ments of the clusters and the grains were taken from indepen- 
dent experiments, see text for details. 

In case of the observed two largely different grain sizes 
this can be expressed as 

-b 3 B 2 s (B,T) (5) 

where B$ = {S z ) /S is the Brillouin- function defined 
under Eq. [3] for the small iron clusters and L(x) — 
cotha; — 1/x is the Langevin- function, i.e. the classical 
limit of Bg for the large Fe grains. The parameters bx, &2 
and &3 represent the relative weights of scattering from 
grain to grain, between a grain and a cluster and from 
cluster to cluster. 

Figure [5] shows magnetoresistance for sample A in the 
superparamagnetic phase and the expected variation cal- 
culated by Eq. |SJ The two characteristic magnetic mo- 
ments (/ig = 500 /iB and S = 17 /ib ) are determined in 
the previous analysis from independent experiments: the 
field and temperature dependence of the bulk magneti- 
zation, and the temperature dependence of the resistiv- 
ity. The relative weights of the various processes were 
used as fitting parameters: b\ = 0.035, 62 = 0.013 and 
63 = 0.0001. The value of bx, i.e the relatively weight for 
the grain— >grain scattering process is larger than &2, even 
though the amplitude of the grain scattering is small£ It 
reflects the large probability of scattering from grain to 
grain due to the large volume fraction of this type of mag- 
netic scatterer. The direct interplay between the clusters 
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is negligible, as expected for a small fraction (less than 
2%). 

In Fig. El the dashed and the dashed-dotted lines rep- 
resents the contribution of the two dominant scattering 
processes. The grain— >grain scattering process is the 
leading term in Eq. up to 4T magnetic field. The 
dashed line corresponds to the square of the magneti- 
zation, and if this would be the only scattering process 
the simple relation^ of AR(B) oc M 2 would hold. Above 
B = 4T, however, the second term of Eq. [5] dominates 
and the observed behavior can be attributed to the ef- 
fect of the small clusters, not seen in the magnetization 
curves. The good description of the measured magneti- 
zation and resistance curves by two characteristic sizes 
is the consequence of the narrow size distributions of the 
grains and the clusters. 

Bimodal distribution of the grain size has already 
been observed in granular systems prepared by co- 
depositionifl* 2 ^ 2 ^, rapid quenching from the melt^S, or 
layered growth 3132 of the constituents. Since for our 
samples a non-saturating magnetoresistance indicates the 
presence of small clusters even in case of 25 nm thick 
continuous Fe layers^ 3 -, we associate the large grains and 
the small culsters to Fe rich grains of the granular lay- 
ers and small Fe clusters trapped inside the Ag layers, 
respectively. Intermixing of the layers can occur dur- 
ing the sample growth even when the heat of mixing is 
positive 2 ^, like in the case of Fe and Ag. On the other 
hand, the tendency for non-equilibrium mixing does not 
seem to depend on the sample preparation method, since 
the magnetoresistance and the magnetic properties of our 
granular multilayers are very similar to those observed in 
co-depositied 2 ^ 3 ^ samples. 

In our case for sample A the magnetic moments of 
the large grains and the small clusters differ more than 
an order of magnitude; /iq =500 /ib, and S w 17/iB- Ac- 
cording to the Mossbauer spectroscopy measurements the 
small clusters contain only a small fraction (below 2 %) 
of the magnetic atoms. This explains why the bulk mag- 
netization is determined by the properties of the large 
grains at all temperatures, while in the transport prop- 



erties the magnetic scattering of small clusters also play 
an important role. 

In case of sample B and C a broader grain size distri- 
bution is indicated by the smeared out ZFC curves and to 
describe the magnetic field dependence of the resistance 
would require further parameters. Interactions between 
the grains are also likely to play a roleii as the Fe concen- 
tration increases. However, the temperature dependence 
of the resistivity measured in 12 T magnetic field is 
quite similar to that observed in sample A (see Fig. 4) 
and the analysis of the magnetic scattering contribution 
(Fig. 8) undoubtedly indicates that there is a significant 
contribution from small clusters in these samples, as well. 

IV. SUMMARY 

In conclusion, we investigated the magnetic scatter- 
ing processes in sequentially evaporated granular Fe-Ag 
films. Unusual magnetotransport features - like sublin- 
ear temperature dependence of the resistivity over a wide 
temperature range both in zero and 12 T magnetic fields 
and large, non saturating GMR - were found experimen- 
tally. The contribution of the magnetic scattering was 
separated and analyzed. The quantitative description 
suggests a granular system with bimodal size distribution 
of the magnetic components: coexisting large grains and 
small clusters. A detailed numerical analysis was given to 
determine characteristic grain- and cluster-moments, and 
their influence on both the scattering processes and on 
the macroscopic magnetization. The analysis reveals that 
scattering on the small clusters plays a dominant role in 
high magnetic fields over a wide temperature range. 
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